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We consider a Taylor-Dean-type flow of an electrically conducting liquid in an annulus between 
two infinitely long perfectly conducting cylinders subject to a generally helical magnetic field. The 
cylinders are electrically connected through a remote, perfectly conducting endcap, which allows 
a radial electric current to pass through the liquid. The radial current interacting with the axial 
component of magnetic field gives rise to the azimuthal electromagnetic force, which destabilizes the 
base flow by making its angular momentum decrease radially outwards. This instability, which we 
refer to as the pseudo-magnetorotational instability (MRI), looks like an MRI although its mech- 
anism is basically centrifugal. In a helical magnetic field, the radial current interacting with the 
azimuthal component of the field gives rise to an axial electromagnetic force, which drives a longi- 
tudinal circulation. First, this circulation advects the Taylor vortices generated by the centrifugal 
instability, which results in a traveling wave as in the helical MRI (HMRI). However, the direction 
of travel of this wave is opposite to that of the true HMRI. Second, at sufficiently strong differential 
rotation, the longitudinal flow becomes hydrodynamically unstable itself. For electrically connected 
cylinders in a helical magnetic field, hydrodynamic instability is possible at any sufficiently strong 
differential rotation. In this case, there is no hydrodynamic stability limit defined in the terms of 
the critical ratio of rotation rates of inner and outer cylinders that would allow one to distinguish 
a hydrodynamic instability from the HMRI. These effects can critically interfere with experimental 
as well as numerical determination of MRI. 



I. INTRODUCTION 

The magnetorotational instability (MRI) can account 
for the formation of stars and entire galaxies in the ac- 
cretion disks. For an object to form, the matter circling 
around it has to slow down by transferring its angular 
momentum outwards. The observed accretion rates sug- 
gest the angular momentum transfer in the astrophysi- 
cal disks to be turbulent while the velocity distribution 
in them seems to be hydrodynamically stable. A pos- 
sible solution to this problem was suggested by Balbus 
and Hawley [l], 0] , who pointed out that a Keplerian ve- 
locity distribution in accretion disk can be destabilized 
by a magnetic field in the process known as the MRI 
[3, Q. This proposition has triggered a number of ex- 
perimental studies trying to reproduce MRI in labora- 
tory 0,0]. The main technical difficulty to such exper- 
iments is the magnetic Reynolds number Rm that is re- 
quired to be ~ 10 at least. For a typical liquid metal 
with the magnetic Prandtl number Pm ~ 10 -5 — 10~ 6 , 
this corresponds to a hydrodynamic Reynolds number 
Re = Rm/Pm - 10 6 - 10 7 0|. Thus, the base flow 
on which the MRI is to be observed may be turbulent 
at such Reynolds numbers independently of MRI as in 
the experiment of Sisan et al. A possible solution 

to this problem was proposed by Hollerbach and Riidi- 
ger @, who suggested that MRI can take place in the 
Taylor-Couette (TC) flow at Re ~ 10 3 when the im- 
posed magnetic field is helical rather than purely axial 
as in the classical case. Theoretical prediction of this 
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new type of helical MRI (HMRI) was soon succeeded by 
a claim of its experimental observation by Stefani et al. 
0, [lfl EH- Subsequently, these experimental observa- 
tions have been questioned by Liu et al. [12j who find 
no such instability in their inviscid theoretical analysis of 
finite length cylinders with insulating endcaps. In a more 
realistic numerical simulation, Liu et al. [13J confirm the 
experimental results, though note that there is no MRI 
at the experimental parameters when ideal TC boundary 
conditions are used. Szklarski [3| showed later that the 
ideal TC requires a slightly different parameters for the 
HMRI to set in. Despite the numerical evidence, Liu et 
al. [3, [H| suspected the observed phenomenon to be a 
transient growth rather than a self-sustained instability. 
This paper shows that the observation of a self-sustained 
instability which looks like an MRI does not necessarily 
mean that the latter is MRI. 

Recently, we found that HMRI can be self-sustained 
and thus experimentally observable in a system of suffi- 
ciently large axial extension because there is not only con- 
vective but also absolute HMRI threshold [if}]. However, 
the comparison with the experimental results [3, [T3, EH 
revealed that HMRI has been observed slightly beyond 
the range of its absolute instability, where it is expected 
to be self-sustained according to the ideal TC flow model. 
This discrepancy with the experimental observations is 
probably due to the deviation of the real base flow from 
the ideal TC flow used in the theoretical analysis. Such a 
deviation, however, poses a major problem for the inter- 
pretation of experimental results, especially for the iden- 
tification of HMRI. Namely, the Rayleigh line defining 
the hydrodynamic stability limit of the ideal TC flow is 
used as a reference point to discriminate between a mag- 
netically modified Taylor vortex flow and HMRI. The 



latter two are hardly distinguishable by the oscillation 
frequency, which varies weakly as the Rayleigh line is 
crossed. The main problem is the hydrodynamic stabil- 
ity limit of the real base flow, i.e., its actual Rayleigh 
line, which may differ from that of the ideal TC flow. 
Therefore the latter cannot be used for the interpretation 
of experimental results. This ambiguity is not resolved 
by the direct numerical simulation of the problem either 
even if there is a perfect agreement with the experiment. 
It is because the notion of MRI is based on the ideal TC 
flow with a fixed hydrodynamic stability limit, which is 
affected neither by the end effects nor by the magnetic 
field. Unfortunately, this is the case neither for experi- 
ments nor for numerical simulations. First, there is an 
Ekman pumping at the endcaps, which can spread up 
to significant but nevertheless limited distances into the 
base flow provided that the latter is hydrodynamically 
stable. The Ekman circulation can be reduced by using 
several independently rotating rings for the endcaps |6| 
or by splitting the latter into two rings of a certain size 
attached to the inner and outer cylinders, respectively 
[T3 |. Another important effect pointed out by Szklarski 
and Riidiger [T3|, which can significantly affect the base 
flow, is related with the Hartmann layers forming at the 
endcaps in axial magnetic field. 

In this paper, we show that there may be additional 
effects in the presence of a magnetic field when well con- 
ducting inner and outer cylinders are electrically con- 
nected through an endcap as in the original PROMISE 
experiment [fj, [1(3, EH. The endcap acting in parallel 
with the Hartmann layer allows a radial current to close 
through the liquid between the cylinders. The interaction 
of radial current with axial magnetic field gives rise to an 
azimuthal electromagnetic force, which reduces the veloc- 
ity difference between the endcap and the liquid above it. 
Depending on the strength of the magnetic field, this elec- 
tromagnetic force can render the profile of azimuthal base 
flow centrifugally unstable. As a result, in axial magnetic 
field, the instability can extend significantly beyond the 
Rayleigh line similarly to the classical MRI. Moreover, 
in helical magnetic field, the interaction of radial current 
with the azimuthal component of magnetic field gives rise 
to an axial electromagnetic force, which drives a longi- 
tudinal flow. First, this longitudinal flow going upwards 
along the inner cylinder, where the azimuthal base flow 
is centrifugally destabilized, advects Taylor vortices that 
results in a traveling wave as in the HMRI. However, the 
direction in which these Taylor vortices are advected is 
opposite to the direction of travel of true HMRI wave. 
Second, for sufficiently large differential rotation, longi- 
tudinal flow may become linearly unstable at any rotation 
rate ratio. 

The paper is organized as follows. In Sec. Owe formu- 
late the problem in the inductionless approximation. The 
base flow for electrically connected cylinders is derived in 
Sec. IIIII Section |IV] introduces the linear stability prob- 
lem. Numerical results for axial and helical magnetic 
fields are presented in Sees. IV Al and IV respectively. 
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Figure 1: (a) An ideal electrically uncoupled system and (b) a 
real system with the inner and outer cylinders electrically con- 
nected via the endcap and inner vessel wall; (1) outer cylinder, 
(2) inner cylinder, (3) inner wall, (4) endcap. 

The paper is concluded with a summary in Sec. IVII 

II. PROBLEM FORMULATION 

Consider an incompressible fluid of kinematic viscosity 
v and electrical conductivity a filling the annulus between 
two long, concentric cylinders with inner radius Ri and 
outer radius R a rotating with angular velocities fij and 
fl . The flow is subject to a generally helical, steady ex- 
ternal magnetic field Bo = B^e^ + B z e z with axial and 
azimuthal components B z = Bq and — (3B^Ri/r in 
cylindrical coordinates (r,<j>,z), where (3 is a dimension- 
less parameter characterizing the geometrical helicity of 
the field. The fluid is supposed to be poorly conduct- 
ing so that the induced magnetic field is negligible with 
respect to the imposed one. This corresponds to the so- 
called inductionless approximation, which holds for the 
HMRI characterized by small magnetic Reynolds number 
Rm = fio<JVQL <C 1, where fj,o is the permeability of vac- 
uum, vq and L are the characteristic velocity and length 
scale [18J . The velocity of fluid flow v is governed by the 
Navier-Stokes equation with electromagnetic body force 

-^ + (v-V)v= — V P + i/V 2 v + -jxB„, (1) 

at p p 

where the induced current follows from Ohm's law for a 
moving medium 

j = a(E + vxB ). (2) 

In addition, we assume that the characteristic time of 
velocity variation is much longer than the magnetic dif- 
fusion time, To ^> T m = po<jL 2 , that leads to the quasi- 
stationary approximation, according to which V x E = 
and E = — V$, where <& is the electrostatic potential. 
Mass and charge conservation require V • v = V • j = 0. 
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III. BASE STATE 

An ideal, axially unbounded system shown in Fig. [0(a) 
admits a translationally invariant base state with purely 
azimuthal velocity distribution Vo(r) = e^voir). Such a 
flow in axial magnetic field induces a radial electric field, 
which gives rise only to the potential difference between 
the inner and outer cylinder but no radial current is in- 
duced because of the charge conservation. Thus, in an 
ideal system, the magnetic field affects the stability of 
the base flow without altering the latter, which is the 
main premise underlying MRI. In reality both cylinders 
may not be completely electrically decoupled from each 
other. For example, in an axially bounded system such a 
coupling may be provided by an electrically conducting 
endcap serving as a closing circuit between the inner and 
outer cylinders. This corresponds to the principal setup 
of the PROMISE experiment illustrated in Fig. O^b). 
The liquid metal in the narrow gap separating the inner 
cylinder (2) from the endcap (4) and the inner wall (3), 
which both form a solid well-conducting vessel together 
with the outer cylinder, provides a sliding contact. This 
allows a radial electric current to pass through the liquid 
between the outer and inner cylinder and then to close 
either directly through the endcap or via the inner wall as 
sketched in Fig. \V[b). Note that this setup is analogous 
to the homopolar generator also known as the Faraday 
disk. 

In the following, an axially uniform radial current 
Jo = e rjo{ r ) is supposed to pass through the liquid and 
close through a remote endcap. Note that this current 
is generated by the differential rotation of cylinders in 
axial magnetic field rather than applied externally as it 
is planned in the so-called Kurchatov MRI experiment 
PJJ, [2fl|- Our main assumption is that the system is suf- 
ficiently extended so that an axially uniform base state 
can develop sufficiently far away from the ends as in the 
classical TC setup. Thus, we neglect any direct effect of 
the endcap on the base flow, which is affected only by an 
axially uniform radial current passing through the liquid. 
The charge conservation yields jo(r) = Jo/r, where Jo is 
a constant that will be determined later by specifying 
the connection between the cylinders. First, the interac- 
tion of radial current with axial magnetic field gives rise 
to the azimuthal electromagnetic force, which affects the 
profile of azimuthal velocity. The latter is governed by 
the <j) component of Eq. |(3j) 
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whose solution can be written as «o(r) — v~a(r) — 
Jo-Bo / (2^ p)vo(r) , where v (r) and Vo{r) are the profiles 
of the classical Couette and electromagnetically driven 
Dean [2l| flows, 
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Linear stability of such a Taylor-Dean (TD) flow in 
purely axial magnetic field has been considered by Szk- 
larski and Riidiger [17l ]. In contrast to us, they regard 
the Dean component of the flow to be independent of the 
Couette one, but ignore that the current driving the for- 
mer is induced by the latter, i.e., the differential rotation 
of the cylinders. 

Second, in a helical magnetic field, radial current in- 
teracting also with the azimuthal component of magnetic 
drives a longitudinal flow wq(z) governed by the z com- 
ponent of Eq. J!} , 
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The solution can be presented as 

where (3% = f3Ri and Wo,i(r) and wo.2(r) are the parts of 
flow driven by the pressure gradient and by electromag- 
netic force 
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The axial pressure gradient Po = dpo/dz, which is con- 
stant for a longitudinally uniform flow, is related to 
the electromagnetically-driven part of the flow by the 
flow rate conservation J„° w (r)rdr = yielding P — 
—2JpBof3iKo, where 
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results in a simple but long analytic expression which 
is skipped here. Eventually, we obtain wp(r) = 
J B f3i/(2pv)w {r), where w Q {r) = K Wo,i(r) - Wo^r) 
depends on the geometry only. In order to determine the 
last unknown quantity Jo, we need to specify how the 
inner and outer cylinders are connected to each other by 
the endcap. In the following, we focus on the experimen- 
tal configuration shown in Fig. QIb), where the outer 
cylinder (1) forms a solid body together with the end- 
cap (4) and inner wall (3), while the inner cylinder (2) is 
separated from the endcap and the inner wall by a rela- 
tively thin gap filled with the liquid metal, which serves 
as a sliding contact. First, we integrate Ohm's law [Eq. 
([2l)] over the liquid gap giving us the radial voltage drop 
between the inner and outer cylinders 
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where T = J R ° [vo(r) + 0f wo(r)/r\ dr represents an- 
other long analytic expression. Second, since there is no 
axial voltage drop along perfectly conducting cylinders, 
the same radial voltage drop can be found alternatively 
by integrating Ohm's law radially from Ri to R a over the 
endcap, which is also assumed to be perfectly conducting, 

*o - * 4 = \b o VL (Rl - Rf) + J K, (4) 

where TZ is a phenomenological parameter introduced to 
account for effective linear resistance of the sliding liquid- 
metal contact between the inner cylinder and the endcap. 
Substituting this into Eq. J3]) we obtain 

which is the last quantity defining the base state. 

Now it remains to estimate the resistance TZ introduced 
in Eq. (J4j) for the setup shown in Fig. |Tfb) that is de- 
scribed in detail in Refs. and As seen, there are 
two parallel paths for the current to connect between the 
inner cylinder (2) and the endcap (4). First, the current 
can connect directly over the vertical gap of w 10 mm 
width between the inner cylinder and the endcap. Sec- 
ond, the current can connect over the annular gap of 
4 mm width between the inner cylinder (2) and the inner 
wall (3) and then pass along the latter towards the endcap 
(4) . Because of a much larger contact area, the effective 
resistance of the second path is obviously much smaller 
than that of the first one, which thus may be neglected 
in this parallel connection. On the other hand, the gap 
width of the second path is by an order of magnitude 
smaller than the 40 mm width of whole liquid layer be- 
tween the inner and outer cylinders. Thus, the resistance 
of the second path may be neglected with respect to that 
of the whole liquid layer, which is connected in series with 
the latter. In the following, we assume TZ = that sup- 
poses a negligible contact resistance between the inner 
and outer cylinders, which appears to be a good approx- 
imation to this PROMISE setup. The limit of TZ -> oo 
corresponds to the classical case of electrically decoupled 
cylinders. Note that in Eq. J5]) TZ stands next to the 
electromagnetic term ~ B 2 implying that even a finite 
TZ may become negligible in sufficiently strong magnetic 
field. In addition, note that the actual PROMISE setup 
is considerably more complex than this simple model. 
In particular, we assume that the sidewalls are perfectly 
conducting with respect to the liquid metal, whereas the 
conductivity of Copper sidewalls is only 13 times higher 
than that of the GalnSn eutectic alloy used in the ex- 
periment. Although our model is relatively rough, it can 
still highlight some principal effects overlooked by more 
elaborate numerical models. 



IV. PERTURBED STATE 

We consider a perturbed state 

{r.'}(""-{2:«}« + { ?,€}('■«' 

where Vi, p%, j 1( and 4>i present small-amplitude pertur- 
bations for which Eqs. HJ and ^ after linearization take 
the form 

<9vi 

— — + (vi • V)v + (v ■ Vjvi 

ot 

1 o 1 

= — Vpi + ^V 2 vi + -j x x B , (6) 

P P 
jj = ct(-V*i+vi xB ). (7) 

In this paper, we focus on axisymmetric perturbations, 
which are typically more unstable than nonaxisymmetric 
ones for TC flow [22], however this is not always the case 
for the conventional TD flow [23]. Analysis of nonaxisym- 
metric perturbations for an electromagnetically driven 
TD flow is outside the scope of this paper. In the ax- 
isymmetric case, the solenoidity constraints are satisfied 
by meridional stream functions for fluid flow and electric 
current as 

v = ve$ + V x (ipej,), j = je^ + V x (he,/,). 

Note that h is the azimuthal component of the induced 
magnetic field, which is used subsequently instead of 
$ for the description of the induced current. Thus, 
we effectively retain the azimuthal component of the 
induction equation to describe meridional components 
of the induced current, while the azimuthal current is 
related explicitly to the radial velocity. In addition, 
for numerical purposes, we introduce also the vorticity 
u) = V x v = uuej, + V x (ve$) as an auxiliary variable. 
The perturbation is sought in the normal-mode form 

{»i,wi,ik,M(r,t) = {v,u,ip,h](r)e^ +ik!1 , 

where 7 is, in general, a complex growth rate and k is 
the axial wave number. Henceforth, we proceed to di- 
mensionless variables by using Ri, Rj^jv, Ri^li : Bq 7 and 
uBoRiVli as the length, time, velocity, magnetic field, 
and current scales, respectively. The nondimensionalized 
governing equations are 

"fv — Dki — Rei/c (^wqv — r _1 (tvq)' tfj + Ha 2 ikh, (8) 
jlu = DkQj — Keik (wquj + r (r~ 1 w' )' ip — 2r~ 1 vo^ 

-Ha 2 ifc (iki> + 2f3r~ 2 hj , (9) 
= D k ijj + u, (10) 
= D k h + ik (v - 2/3r~ 2 ?/^ , (11) 

where D^f = r _1 (rf) — (r~ 2 + k 2 )f and the prime 
stands for d/dr; Re = R 2 Q,i/v and Ha = RiBo^J a / (pv) 
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Figure 2: Angular momentum profiles for (a) = 0, (b) /i = 0.25, (c) /i = 0.5 and various Hartmann numbers for cylinders 
with A = 2 in axial magnetic field (/3 = 0) with no contact resistance (TZ = 0). 



are the Reynolds and Hartmann numbers, respectively. 
Boundary conditions for the flow perturbation and the 
electric stream function on the perfectly conducting inner 
and outer cylinders at r = 1 and r = A, respectively, are 
v = $ = = (rh)' = 0. 

The governing Eqs. (f8|- fTl"j) for perturbation ampli- 
tudes were solved in the same way as in Refs. [3] and [IB] 
by using a spectral collocation method on a Chebyshev- 
Lobatto grid with a typical number of internal points 
N = 32, which ensured the accuracy of about five digits. 

The dimensionless azimuthal and axial velocity com- 
ponents of the base flow 

v {r) = v (r) + -Ea 2 J v Q (r), (12) 
w (r) = iHa 2 J /Mo(r), (13) 

follow straightforwardly from the corresponding dimen- 
sional counterparts when iij and f2j are replaced by 1, 
R a and Q by A = R Q /Ri and \i = o /0j, respectively, 



and /3i by 0; the dimensionless counterpart of Jo is 

x rH + ln\+ ^Ha 2 Z^ , (14) 

where TZ = TZa and 1 = 1jR\ are the dimensionless 
counterparts of TZ and X, respectively. As seen from 
Eqs. (fl2jl- lfl4| . for Ha > 1 velocity profiles tend to 
asymptotic solutions which, as noted above, depend nei- 
ther on the contact resistance TZ nor on the magnetic field 
strength 

v (r) w v (r) + J v (r), (15) 
w (r) w J /3w (r), (16) 

where J = (/i - 1) + i) /X 

V. NUMERICAL RESULTS 
A. Axial magnetic field 

We start with an axial magnetic field (f3 = 0), for which 
the base flow is purely azimuthal. The profiles of angular 



6 




0.1 0.2 0.3 ,, 0.4 0.5 0.6 0.7 0.1 0.2 0.3 ,, 0.4 0.5 0.6 0.7 



Figure 3: (a) Critical Reynolds number Re c and (b) wave number k c versus fi, at various Hartmann numbers for cylinders with 
A = 2 and no contact resistance (1Z = 0) in purely axial magnetic field. 



momentum rvo(r) are shown in Fig. [2] for several cylin- 
der rotation rate ratios fi and various Hartmann num- 
bers. For fx = 0, shown in Fig. OJa), which corresponds 
only to the inner cylinder rotating, the profile without 
the magnetic field is centrifugally unstable with the an- 
gular momentum decreasing radially outward. In this 
case, the magnetic field slows down the overall rotation 
rate of the liquid making the angular momentum decrease 
faster at the inner cylinder that may result even in the 
reversal of the sense of liquid rotation at the outer cylin- 
der when the magnetic field is sufficiently strong. This 
effect is due to the magnetic field trying to eliminate the 
differential rotation between the liquid and the endcap, 
which is attached to the outer cylinder and thus rotates 
with a lower angular velocity than the liquid above it as 
long as < 1. A similar effect can also be observed in 
Fig. [2jb) for /j, = 0.25, which without the magnetic field 
corresponds to a marginally stable state with a constant 
angular momentum distribution. In this case, the mag- 
netic field again retards the liquid rotation so rending the 
distribution of angular momentum centrifugally unstable 
at the inner cylinder and stable at the outer one. For 
/i = 0.5 shown in Fig. Etc), the profile without magnetic 
field is centrifugally stable with the angular momentum 
increasing radially outward. However, a strong enough 
magnetic field changes the distribution of the angular 
momentum at the inner cylinder from radially increas- 
ing to decreasing one so rending the profile centrifugally 
unstable. 

This is confirmed by the critical Reynolds number plot- 
ted against [i for various Hartmann numbers in Fig. [3]^ a) 
with the corresponding critical wave numbers shown in 
Fig. [3l[b). As seen in Fig. [HJa), without magnetic field 
(Ha = 0), the critical Reynolds number tends to infinity 
as /i approaches the Rayleigh line ji c = A~ 2 = 0.25 de- 
fined by d (r 2 f2) jdr = 0, at which the profile of angular 
momentum becomes centrifugally stable. As the mag- 
netic field is increased, the instability starts to extend 
beyond the Rayleigh line reaching [i ~ 0.65 at sufficiently 



large Hartmann numbers. Although this extension of 
the instability beyond the Rayleigh line may look like an 
MRI, it has a principally different physical mechanism. 
Namely, in the MRI, the magnetic field destabilizes the 
flow without altering it, whereas here the magnetic field 
does alter the base flow by rendering it centrifugally un- 
stable as discussed above. Moreover, the standard MRI 
in axial magnetic field is not captured by the induction- 
less approximation (Pm = 0) used here [24]. Thus, in 
axial magnetic field, this centrifugal instability occurring 
beyond the Rayleigh line can easily be distinguished from 
the true MRI. 



B. Helical magnetic field 

As seen in Fig. IH in a helical magnetic field, the base 
flow besides the azimuthal component has also an ax- 
ial one, which is driven by the interaction of radial cur- 
rent with the azimuthal component of magnetic field. In 
the configuration with the endcap attached to a slower- 
rotating outer wall, the induced electric current is flow- 
ing radially outward, as discussed above and, thus, the 
resulting axial electromagnetic force is directed upward. 
Because both the current and azimuthal magnetic field 
decrease radially outward as ~ 1/r, the resulting elec- 
tromagnetic force is stronger at the inner wall, where it 
drives the liquid upward as seen in Figs. IU(b) , [4]Jd) , and 
IDJf). Since the annular gap is supposed to be closed at 
both ends, the constant axial pressure gradient arising 
in the response to the electromagnetic force drives a re- 
turn flow along the outer cylinder, which compensates 
for the upward one along the inner cylinder. This ax- 
ial flow in the azimuthal magnetic field, in turn, induces 
an additional electrostatic potential, which contributes 
to that induced by the azimuthal flow in the axial field 
as described by Eq. ([3]). The total potential difference 
induced by the flow between the inner and outer cylin- 
ders balances that induced by the rotation of bottom in 
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Figure 4: [(a), (c), and (e)] Radial profiles of angular momentum and [(b), (d), and (f)] axial velocity for [(a) and (b)] /i = , 
[(c) and (b)] /i = 0.25, [(e) and (f)] \i = 0.5 and various Hartmann numbers for cylinders with A = 2 and no contact resistance 
(TZ — 0) in helical magnetic field with p = 6. 



the axial magnetic field, which is given by Eq. . The 
potential balance determines the magnitude of the in- 
duced radial current defined by Eq. |(5]), which, in turn, 
interacts with the magnetic field and disturbs the flow. 
Thus, the perturbation of the azimuthal flow is weaker 
in helical magnetic field than it is in a purely axial one 
because a part of the potential difference is compensated 
by the axial flow [see Figs. [2 (Ha), gfc), andHJe)]. 

The instability characteristics in a helical magnetic 
field plotted in Fig. [5] differ considerably from those in 



axial magnetic field shown in Fig. [31 In contrast to the 
axial magnetic field, now the most unstable mode of in- 
stability is oscillatory, i.e., a traveling wave as for the 
HMRI. However, it is important to note that the phase 
velocity of this wave, which is determined by the sign 
of the frequency shown in Fig. (SJc), is directed upward 
oppositely to that of true HMRI. The reversed phase ve- 
locity is due to the longitudinal flow, which is absent for 
the ideal HMRI with electrically decoupled cylinders. As 
seen in Figs. [UJa),[Hc), and[D[e), the radial current in- 
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Figure 5: (a) Critical Reynolds number Re c , (b) wave number k c and (c) the rescaled frequency u c /Re c versus fi at various 
Hartmann numbers for cylinders with A = 2 and no contact resistance (H = 0) in helical magnetic field with (3 = 6. 



teracting with the axial component of the magnetic field 
causes the angular momentum to decrease radially out- 
ward at the inner cylinder that renders the flow centrifu- 
gally unstable. Furthermore, the Taylor vortices arising 
at the inner cylinder are advected by the longitudinal flow 
upward. The advection in this case obviously dominates 
over the direct electromagnetic effect of helical magnetic 
field, which would drive the true HMRI wave in the op- 
posite direction. 

Moreover, in a helical magnetic field in contrast to 
purely axial one, the instability is seen to extend much 
farther beyond the Rayleigh line up to the limit of solid- 
body rotation defined by [i = 1 and even beyond it, which 
is not considered here. The instability in helical magnetic 
field differs significantly from that in purely axial field. 
As seen in Fig. \E[ for Ha = 1, shortly after the Rayleigh 
line, the most unstable mode switches from the initial 
Taylor vortices branch to another one, which is obviously 
associated with the axial flow. For larger Hartmann num- 
bers, this transition proceeds smoothly with the critical 
wave number developing a maximum at certain fi when 
Ha < 10. Beyond the Rayleigh line, the critical Reynolds 
number first decreases with the Hartmann number up to 
Ha « 10 and then starts to grow for larger Ha again. 
For larger //, the most unstable mode jumps to another 



branch with a considerably smaller critical wave number 
and positive frequency, which corresponds to the oppo- 
site direction of the phase velocity now coinciding with 
that of the true HMRI. Note that such jumps of the crit- 
ical mode are characteristic also for the conventional TD 

flow dHHi]. 



VI. SUMMARY AND CONCLUSIONS 

We have considered linear stability of a TD-type flow of 
an electrically conducting liquid in the annulus between 
two infinitely long perfectly conducting and differentially 
rotating cylinders in the presence of a generally helical 
magnetic field. The cylinders were supposed to be elec- 
trically connected through a remote endcap. We showed 
that this electrical connection can render the base flow 
hydrodynamically unstable. First, the azimuthal base 
flow in an axial magnetic field gives rise to a radial emf. 
If the cylinders are electrically decoupled, no current can 
close between them, and, consequently, the emf results in 
the radial charge redistribution, which gives rise to the 
electrostatic potential whose gradient compensates the 
original emf. If there is no current, there is no electro- 
magnetic force and no effect of the magnetic field on the 
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base flow either. This corresponds to the ideal TC flow, 
which is used as a reference for the definition of MRI, 
where the magnetic field is expected to destabilize the 
base flow by affecting only its disturbances but not the 
base flow itself. 

This is no longer the case when the cylinders are elec- 
trically connected, and a radial current can close between 
them. The interaction of radial current with the axial 
component of the magnetic field gives rise to the az- 
imuthal electromagnetic force, which tries eliminate the 
velocity difference between the endcap and the liquid 
above it. Depending on the strength of magnetic field 
and the effective contact resistance between the inner and 
outer cylinder, this electromagnetic force can modify the 
profile of azimuthal base flow so that it becomes cen- 
trifugally unstable. As a result, the magnetic field makes 
the instability extend significantly beyond its apparent 
Rayleigh line so resembling MRI in the case of an un- 
perturbed TC flow. Furthermore, in a helical magnetic 
field, the interaction of radial current with the azimuthal 
component of magnetic field gives rise to an axial electro- 
magnetic force, which drives a longitudinal flow. First, 
this longitudinal flow going upward along the inner cylin- 
der, where the azimuthal base flow is destabilized by the 
magnetic field, advects Taylor vortices, so giving rise to a 
traveling wave as in helical MRI. However, the direction 
of the most unstable traveling wave of this centrifugal 
instability is opposite to that of the true MRI. Second, 
for sufficiently large differential rotation, the longitudi- 
nal flow becomes hydrodynamically unstable itself. For 
electrically connected cylinders in helical magnetic field, 
hydrodynamic instability is possible at any sufficiently 
large differential rotation. In this case, there is no pure 
hydrodynamic stability limit defined in the terms of the 
critical ratio of rotation rates of inner and outer cylinders 
that would allow one to discriminate between magneti- 
cally modified hydrodynamic instability and the HMRI. 

From the experimental point of view, a crucial test 
for the pseudo-MRI would be the extension of the Tay- 
lor vortex flow beyond the Rayleigh line in purely axial 
magnetic field at Rm < 1. The PROMISE experiment re- 
ports only one such apparently successful test in which, 



however, the time-averaged flow and thus stationary Tay- 
lor vortices, if any, are removed [§]. Traveling wave ap- 
pears as soon as the azimuthal component of the field is 
switched on. As to the helical magnetic field, the exper- 
iment seems to find the right direction of the phase 
velocity in agreement with the ideal HMRI model rather 
than that of the pseudo-MRI considered in this paper. 
But this does not necessarily mean that the real base 
flow in the experiment is any closer to the ideal TC one. 
Note that the nonaxisymmetric m = 1 instability mode 
unexpectedly observed in the PROMISE experiment is 
characteristic for certain regimes of the conventional TD 

flow H3. 

Although the current circulation through the liquid 
metal has been eliminated in a modified PROMISE 
experiment[27] by insulating the inner cylinder, the base 
flow still remains strongly affected by the Ekman pump- 
ing due to the endcaps which makes it more complex 
than the one used in this study. In the new PROMISE- 
2 setup [281 . the Ekman pumping has been signifi- 
cantly reduced by using split rings for the endcap, which 
is insulating now, and thus prevents the current circula- 
tion through it. Although the instabilities appear much 
sharper in the new setup than in the previous one, the 
actual hydrodynamic stability limit, if any, of the base 
flow and so the nature of the observed instabilities is still 
unclear. In particular, as shown by Szklarski and Riidi- 
ger the base flow may significantly be affected by 
the magnetic field also when the endcaps are insulating 
provided that Ha > 10. 

In conclusion, it is not appropriate to use the Rayleigh 
line of the ideal TC flow as a criterion to determine the 
MRI in a significantly different flow. More elaborate nu- 
merical analysis may be necessary for this purpose. 
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